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Abstract 

We propose a scalable scheme to implement quantum computation in graphene nanoribbon. It 
is shown that electron or hole can be naturally localized in each zigzag region for a graphene 
nanoribbon with a sequence of Z-shaped structure without exploiting any confined gate. An one- 
dimensional graphene quantum dots chain is formed in such graphene nanoribbon, where electron 
or hole spin can be encoded as qubits. The coupling interaction between neighboring graphene 
quantum dots is found to be always-on Heisenberg type. Applying the bang-bang control strategy 
and decoherence free subspaces encoding method, universal quantum computation is argued to be 
realizable with the present techniques. 
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Electron spin is one of the leading candidates for the realization of a practical solid 
qubit [lj . The coherent manipulation of electron spins in GaAs quantum dots has been effi- 
ciently realized [2,3]. However due to the interaction with the environment, the decoherence 
time is often in nanoseconds scale in GaAs quantum dots . Even by applying the 

complex technique to prepare nuclear state, the dephasing time for spin qubits is just about 
Ifis |6]. The decoherence is one of the most challenges in the way to quantum computer 
in GaAs quantum dots. Due to the weak spin-orbit coupling and hyperfine interactions 
in carbon, graphene is argued to be an excellent candidate for quantum computation 
However, due to the special band structure of graphene |8|], the low-energy quasiparticles in 
graphene behave as Dirac fermions, and the Klein tunneling and Chiral effect lead to the 
fact that it is non-trivial to form good quantum dot (localized electron states) in graphene. 
It has been shown that the massless Dirac fermions in graphene can ben confined by using 
suitable transverse states in graphene nanoribbons (GNR) , by combining single and 
bilayer regions of graphene 10, III] or by using inhomogeneous magnetic fields 



121. Re- 



cently, there was an experiment report that GNR with well defined zigzag or armchair edge 
structures can be chemically produced fl3l ]. It has also been discovered that localized states 
exist in the zigzag region in Z-shaped GNR 
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Here we present a scalable quantum computation scheme based on Z-shaped GNR quan- 
tum dot system without exploiting any confined gates. The localized particle can be chosen 
to be electron or hole by adjusting the back gate even in the room temperature. The qubit 
is encoded on the electron (hole) spin states localized in the zigzag region of the GNR with a 
sequence of Z-shaped structure. The interaction between qubits is determined by the GNR 
geometrical structure and found to be Heisenberg form. By exploiting bang-bang (BB) con- 
trol strategy and decoherence free subspaces (DFSs) encoding method, universal quantum 
gates are shown to be realizable in this system with the present techniques. 

Based on the 7r orbital tight-binding model, the local density of state (LDOS) and the 
band structure of the zigzag region in a GNR with a sequence of Z-shaped structure can be 
obtained by the direct diagonalization of the single particle Hamiltonian Hq = Ylij r yN)(il> 
where the hopping matrix element = — r if the orbits \i) and \j) are nearest neighboring 
on the honeycomb lattice, otherwise = Q • From the calculated band structure, 
we can see that there are several localized states with electron-hole symmetry around the 
zero energy point as shown in Fig. [TK Thus we can choose to get one localized electron or 
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hole in the zigzag region by adjusting the Fermi level through the back gate. The electron 
ground state energy and the energy gap between the ground state and the first excitation 
state are very sensitive to the size of the zigzag region, as shown in Fig[2J It has been known 
that the width of the armchair GNR (N unit cells) decides whether the system is metallic 

n n 

or semiconducting [15|, |16|. li N = 3m — 1 (m being an integer), the armchair GNR is 
metallic, otherwise it is semiconducting. In addition, for the present Z-shaped structure the 
boundaries along the ribbon of armchair region is unsymmetrical when N is even. Actually, 
in our calculation we find there is no confined state in the zigzag region of Z-shaped GNR 
when N = 3m — 1 or iV = 2m as shown in Fig. [2j On the other hand, when N is 7 and 
the length of the zigzag region L (unit cells) (see Fig. [Tb) is 3, 4, 5, 6, 7 and N = 9, L = 3, 4 
both the ground level and energy gap are above O.leV. Thus we can confine electron (hole) 
to form quantum dot even in the room temperature. 

Fig. [lb shows the spatial distribution of local density of ground state for a GNR with two 
Z-shaped structure in series. Each zigzag region (quantum dot) confines one electron and 
the quantum dots are coupled by the exchange coupling J\. We can obtain J\ by calculating 
the exchange integral J\ = J (p\ (f[ ) ip \ (r~2 ) r <Pi C 7 ^) ^2 (rt) df[ df?. , where (pi(r) and ^(r) 
are the wavefunction of neighboring graphene quantum dots. We can also calculate the next 
nearest neighboring exchange coupling J 2 by the same method. Obviously, the exchange 
coupling Ji, J2 are determined by the geometrical structure of the nanoribbon. For each 
N and L, J\ and J\j J 2 depend on the number of unit cells (D) between two neighboring 
qubits. By numerical calculations, J\, J\j J% are obtained with different N, L, D, as shown 
in Fig. [31 For iV = 7, L = 6, D = 18, J x = 8/ieV, J x j ' J 2 = 10 5 , we can safely neglect this 
non-nearest neighboring qubits coupling. For clarity, in the following discussion we focus on 
the atomic structure with N = 7, L = 6, D — 18. 

To carry out quantum computation, we have to form the logical qubit and realize universal 
quantum gate. It has been shown that single qubit rotations combined with two-qubit 
operations can be used to create basic quantum gates [17]. The spin of the localized electron 
or hole can be used as the physical qubit and the GNR with a sequence of Z-shaped structure 
forms an one-dimensional qubit chain as shown in Fig. [TH. The neighboring qubits in 
this chain have an always-on Heisenberg interaction H = J\S\ ■ Here Si and S2 are 
the spin operator of the neighboring localized electron (hole). It has been known that 
BB control strategy and DFSs encoding method do not require directly controlling the 
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FIG. 1: (Color online) Schematic of the proposed architecture of GNR for quantum computation, 
(a) The Z-shaped GNR quantum dot can localize one electron (left figure) or hole (right figure) in 
the zigzag region by adjusting the Fermi level through the back gate, (b) Local density of states of 
GNR with two Z-shaped structure in series, (c) A GNR with two Z-shaped structure in series, each 
zigzag region confines one electron. The physical qubit is encoded into the spin of the confined 
electron, (d) Special encoding method to eliminate the interaction between logical qubits. Physical 
qubits 1 and 2 form logical qubit L\\ physical qubits 3 and 4 form logical qubit L<i\ physical qubits 
5 and 6 form logical qubit L3. The G z , G x , G y are the BB operation sets of L\, L2 and L3 
respectively. 

interaction between qubits Q] • The quantum information in qubits can be protected 
from decoherence induced by the environment and undesired disturbance induced by the 
inherent qubit-qubit interaction with these strategies. 

For a sequence of Z-shaped structure GNR with N = 7, L = 6, D = 18, the Hamiltonian 
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FIG. 2: (Color online) The ground energy level (Eq) and energy gap (E\ — Eq) between the ground 
state and the first excitation state of the Z-shaped GNR quantum dot device with different width 
of nanoribbon (N unit cells) and different length of quantum dot region (L unit cells). 
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FIG. 3: (Color online) (a) The coupling energy J\ of two nearest neighboring qubits dependence on 
the number of unit cells between two qubits is presented for different size of quantum dot region, 
(b) The ratio of J\ to the next nearest neighboring exchange energy J2 depend on the distance D 
of two neighboring qubits for different size of quantum dot region. 
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of the system can be expressed as 

where (T^f ,z are the spin Pauli operators of the localized electron (hole) in the quantum dots, 
% and j represent two neighboring dots. Here we have neglected the interaction between 
non-neighboring dots, which has been shown to be 5 orders smaller than the neighboring 
interaction. 

To avoid the spin qubits to entangle with the environment, we can apply a BB operation 
U z = exp(—ia z 7r/2) to each quantum dot region. Such rotation operations can be realized if 
a pulsed magnetic field could be applied exclusively To counteract phase decoherence, 



20j. For a simply DFSs encoding, two physical qubits can encode 



we can use DFSs encoding 
a logical qubit: 

|0)l = ITi| 2 ),|i)l = |IiT2). (2) 

As shown in Fig. [lb, we use localized electron in the two neighboring zigzag regions to form 
a logical qubit. 

In order to protect quantum information in the logical qubits, we must decouple the 
always-on Heisenberg interaction between two physical qubits within a logical qubits and 
interaction between two neighboring logical qubits. A nonsynchronous BB pulse operations 
and a special encoding method can be exploited to eliminate these interactions [3] . Here we 
propose an architecture in which the one-dimensional GNR chain form a periodic structure 
L1L2L3L1L2L3 ■ ■ ■ with three logical qubits as a unit, as shown in Fig. [Tbl. L\ represents a 
logical qubit encoded as Eq.(j2]). L 2 is a logical qubit encoded as 

|0)l2 = ~(| t) 3 + l l) 3 )(l T>4 — 1 1)4), (3) 

|1)m = ^(| t>3 — I i> 3 )(| t)4 + | 1)4). (4) 
And L 3 is a logical qubit encoded as 

|0>LS=^(|T>B + <U>5)(|T)6-<U> 6 ), (5) 
|l>L3=i(lt>5-<U>6)(|t>6 + «U>6). (6) 

With this periodic architecture, we have to apply nonsynchronous BB pluse operations 
respectively to L 1; L 2 , L 3 from the operation set G z = {I,U Z , R z }, G x = {I,U X , R x }, 



G y = {I, Uy, R y }, where U z = —o\ ® erf, R z = —ill <8> erf, U x = —of ® erf, R x = ® af , 
Uy = —o\ ® erf, and R y = ® cr|. Then we obtain a quantum computation system with 
entirely decoupled logical qubits. 

Now we show how to carry out universal quantum gates of the logical qubits defined 
above. Logical operations X and Z can generate all SU(2) transformations of logical qubit. 
For logical qubit L\ y X = |(crf ®a 2 + <rf ®cr 2 ), Z = \{o\ — erf). X can be easily achieve d by 



recoupling qubits 1 and 2 by adjusting the BB pulses of both qubits to be synchronous [19]. 
The operation time can be obtained by J At = hir/A , for iV = 7, L = 6 , D = 18, 
At = 0.2ns. Z can be achieved by directly varying the Zeeman splitting on the two physical 
qubits individually analogous to single-qubit operations in the Loss-DiVincenzo quantum 
computer [l|. The operation time of this Z gate can be about Ins when 20mT magnetic 
field could be pulsed exclusively onto each quantum dot region. The fidelity of the X gate can 
be effected by the fluctuation or inhomogeneity in the exchange coupling Ji between different 
dots. The charge noise in the back gate can cause the fluctuation of J%. The main sources of 
the J\ inhomogeneity between different dots are the disorder, irregular edges and defect of the 
GNR. Short-range disorder scarcely affects the LDOS of the ground state. The irregular edge 



effect and long-range disorder change the LDOS but do not destroy the confined states [14]. 
To get high fidelity for X operation, we should avoid the long-range disorder and irregular 
edge. Actually, if we know the coupling Ji between different dots exactly, inhomogeneity 
of Ji can not effect the fidelity of X gate when corresponding inhomogeneous operation 
times are used. In addition, we find the effect of the J\ inhomogeneity or fluctuation to the 
fidelity of the X gate is small as shown in Fig. HI Because the nuclear field would change the 
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evolution of the spin states, the fidelity of the Z gate is dominated by the nuclear field 
The fidelity of the Z gate can be very high due to small nuclear field in graphene system. 
Similarly, high fidelity operation X and Z can be also realized for logical qubits L 2 and L 3 . 
We can construct CNOT gate between two neighboring logical qubits, for example L\ 



and L 2 , by W gate W = |0}(0| ® I + |1}(1| ® e i29Z = e i9Z ® z and Hadamard operation [21] 
By performing Hadamard transformation 
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to the two physical qubits of the second logical qubit L 2 and changing the BB control pulse 
to be the same with L 1; we can recouple the two neighboring logical qubits and implement 
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FIG. 4: The fidelity of X gate against the fluctuation or inhomogeneity of the exchange coupling 
energy J\. 



W gate of logical qubits of L\ and L2. For the present graphene quantum dots chain with 
N = 7, L = 6, D = 18, the total operation time of a CNOT gate can be implemented in 
about Ins with an oscillating magnetic field of lOOmT to achieve the Hadamard operation. 
Similar to the above discussion for X and Z operation, we can find that the fluctuation or 
inhomogeneity of J\ and the nuclear field have trivial effect to the fidelity of the CNOT gate 
in the present protocol. 

The major decoherence sources of spin qubits in solid state system have been identified 
as the spin-orbit interaction and hyperfine interaction. The weak spin-orbit coupling have 
been predicted in carbon material due to the low atomic weight 22j. Since the primary 
component of natural carbon is the zero spin isotope 12 C, the very long coherence time 
given by hyperfine coupling has been theoretically argued [7|. Assuming the abundance of 
13 C is about 1% as in the nature carbon material, the decoherence time has been predicted 
to be more than 10/is in the graphene quantum dot H]- This decoherence time is 4 orders 
longer than the gates operation time of the present protocol. In addition, the decoherence 
time can be much longer if the percentage of 13 C is decreased by isotopic purification. 

In this paper we have presented a scalable scheme of quantum computation based on 
GNR with a sequence of Z-shaped structure. No confined gates is needed to localize the 
particle, which can be chosen to be electron or hole by adjusting back gate. The qubit 
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is encoded in electron or hole spin states, which is naturally localized in the zigzag region 
of GNR even in room temperature. The neighboring qubits are found to have an always- 
on Heisenberg interaction and the dynamical decoupling techniques with DFSs is exploited 
to achieve universal quantum computation in this system. Due to recent achievement in 
production of graphene nanoribbon, this proposal may be implementable within the present 
techniques. 
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